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Abstract. For the Laplacian A defined on a p.c.f. self-similar fractal, we give an explicit 
formula for the resolvent kernel of the Laplacian with Dirichlet boundary conditions, and 
also with Neumann boundary conditions. That is, we construct a symmetiic function G* '' 
which solves {AI - A)"'/(a:) = j G'-^^{x,y)f(y) dfi(y). The method is similar to Kigami's 
construction of the Green kernel in [KigOl, §3.5] and is expressed as a sum of scaled and 
"translated" copies of a ceitain function i^' " which may be considered as a fundamental 
solution of the resolvent equation. Examples of the explicit resolvent kernel formula are 
given for the unit interval, standard Sierpinski gasket, and the level-3 Sierpinski gasket 
SG3. 



Contents 

1 . Introduction 1 

2. The resolvent kernel for the unit interval 5 

3. The Dirichlet resolvent kernel for p.c.f. self-similar fractals 8 

4. The Neumann resolvent kernel for p.c.f. self-similar fractals 15 

5. Example: the Sierpinski gasket SG 16 

6. Example: SGj, a variant of the Sierpinski gasket 19 
References 26 



1. Introduction 

A theory of analysis on certain self-similar fractals is developed around the Laplace 
operator A in [KigOl]. In this paper, we consider the resolvent function {AI - A) ' and 
obtain a kernel for this function when the Laplacian is taken to have Dirichlet or Neumann 
boundary conditions. That is, we construct a symmetric function C " which weakly solves 
{AI - A)G'"(Jc,3') = S{x,y), meaning that 

J G"\x,y)f{y)dfi{y) - (^I - Ar'f{x). (1.1) 
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For the case A — 0, this is just the Green function for A. Consequently, it is not surprising 
that our construction is quite analogous to that of the Green function as carried out in 
[KigOl, §3.5]; see also [Str06, §2.6] for the case of the Sierpinski gasket (and the unit 
interval) worked out in detail, and [Kig03]. 

We present our main results in § 1 .2, just after the introduction of the necessary technical 
terms in §1.1. It is the authors' hopes that the resolvent kernel will provide an alternate 
route to obtaining heat kernel estimates (see [FHK94, HK99]) in this setting, as well as 
other information about spectral operators of the form 



in the same manner as used by Seeley [See67, See69] for the Euclidean situation. Some 
initial results in this direction will appear in [Rog08]. 

To explain the method of construction for the resolvent kernel, we carry out the pro- 
cedure in the case of the unit interval in §2; we believe this particular method has not 
previously appeared in the literature. In §3, we show how the construction may be general- 
ized to any post-critically finite self-similar fractal. In §5, we give the explicit formulas for 
the Sierpinski gasket and in §6 we give the explicit formulas for a variant of the Sierpinski 
gasket which we call SG3. 

Acknowledgements. The authors are indebted to the referee for many keen observations 
and a scrupulously detailed report. 

1.1. Background, notation, and fundamentals. We work in the context of post-critically 
finite (p.c.f.) self-similar fractals. The full and precise definition may be found in [KigOl, 
Def. 1.3.13], but for the present context it suffices to think of such objects as fractals which 
may be approximated by a sequence of graphs, via an iterated function system (IPS). A 
more general setting is possible; cf. [Kig03]. We now make this more precise. 

Definition 1.1. Let {Fi,F2, . . .,Fj}he a collection of Lipschitz continuous functions on 
W' with < Lip{Fj) < 1 for each j. Let X denote the attractor of this IPS, that is, X 
is the unique nonempty and compact fixed point of the set mapping F{A) :- Ijj^j Fj{A). 
The set X is frequently also called the self-similar set associated to this IPS; existence and 
uniqueness of X was shown in [HutSl]. 

Prom the IPS introduced in the previous definition, we now build a sequence of graphs 
which approximates X in a suitable sense. 

Definition 1.2. We use o) - (jD\U2 ■ ■ ■ com to denote a word of length \to\ - m on the symbol 
alphabet {1, 2, . . . , 7). This notation is used to denote a composition of the mappings Fj 
via F^ - Foj^ o Fij, o . . . o . Similarly, K^j - F^{X) refers to a certain m-level cell. The 
collection of all finite words is denoted W» := UmU' 2, . . . , J]"' . 

Definition 1.3. Each map F j of the IPS defining X has a unique fixed point Xj. The 
boundary of X is the largest subset Vq Q {x\, . . . , xn) satisfying 

F^{X) n F^,{X) c F^(yo) n F^,{Vo), for any oj * u' with \u\ = 

The p.c.f. condition mentioned above means that cells Fj{X) intersect only at points of 
FjiVo). The boundary of an m-cell is dK^ :— Fi^{Vo). 

Let Go be the complete graph on Vo, and inductively define G,„ := F{Gm-\)- Also, we 
use the notation x ~,„ y to indicate that x and y are m-level neighbours, i.e., that there is 
an edge in Gm with endpoints x and y. We use Vm - F'"{Vo) to denote the vertices of G,„, 




(1.2) 
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and y* := {j„, V,,,- The fractal X is the closure of V, with respect to either the Euclidean or 
resistance metric. A discussion of the resistance metric may be found in [Str06, § 1 .6] or 
[KigOl, §2.3]. 

Now we are able to make precise the sense in which X is the limit of graphs: one 
may compute the Laplacian (and other analytic objects, including graph energy, resistance 
distance, etc.) for functions m : X — > R by computing it on Gi„ and taking the limit as 
m — > oo. 

Definition 1.4. We assume the existence of a self-similar (Dirichlet) energy form S on X. 
That is, for functions m : X — > R, one has 

J 

£(M) = 2rj'£(MoF,), (1.3) 

for some choice of renormalization factors r\, . . . ,rj e (0, 1) depending on the IPS. This 
quadratic form is obtained from the approximating graphs as the appropriately renormal- 
ized limit of £g„,(«) '■- ^c,„{u, u), where the w-level bilinear form is defined 

£g,„(m, v) := i Yj ^■^>'(«(^) - «0'))(Hx) - v(y)). (1.4) 

The constant C jy - c""l} refers to the conductance of the edge in G,„ connecting x to y (with 
c„ = if there is no such edge). The dependence of c^"l' on m is typically suppressed, as 
X ~y for at most one value of m on p.c.f. fractals. The domain of & is 

domfi := [u : X -> R -£(m) < oo). 

Definition 1.5. We also assume the existence of a self-similar measure /i 

J 

7=1 

with weights fij satisfying < //j < 1 and Y^jf^j - !> ™d normalized so that /i{X) = 1. 
With the notation of Definition 1.2, the measure of the m-cell K^^ is denoted by fiiK^) - 
A'w '■— jUwijUfj, ■ ■ -fJ^ai,,,- The standard measure refers to the case fij — j, for each j. 

Remark 1 .6. The renormalization factor ry should be confused neither with the contraction 
factors Lip{F j) of the maps of the IPS, nor the weights /iy of the self-similar measure yU. 
The values of these constants are completely independent. 

Also, it should be noted that the existence of a self-similar energy asserted in Def- 
inition 1.4 is a strong assumption. While the the self-similar measures of Definition 1.5 
always exist [Hut8 1], the existence of the self-similar energy is a much more delicate ques- 
tion; cf. [Sab97]. 

Definition 1.7. The Laplacian is defined weakly in terms of the energy form. Por u e 
dom£ and / continuous, one says u € dom A with Am = / iff 

&{u,v) - - I fvdp, for all V e domo £, (1-6) 
Jx 

where domo£ is the set of functions in dom£ which vanish on dX - Vo. Note that the 
Laplacian depends on the choice of measure p. 
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More generally, if (1.6) holds with / e L^{d^), then one says u e dom^2 A; and if 

&{u,v) - - I vt/yu, for all V e domo fi, (1-7) 
Jx 

for a finite signed measure yu with no atoms, then one says u e dom^ A. 

It follows from (1.3), (1.5) and Definition 1.7 that A satisfies the scaling identity 

A(m o Fj) = rjHj(hu) o Fj, (1.8) 
and pointwise formula given by the uniform limit 

^u{x) = Um ( f /if ^c///] A„,m(jc), forx e V, \ Vo, (1.9) 

where /if is a piecewise harmonic spline satisfying h^r\y) - ^xy foiy e V,„, and 

A,„u(x) - ^ c^y{u(y) - u(x)), for x € V,„. (1-10) 

V~.f 

Roughly speaking, /if is a "tent" function with peak at x which vanishes outside the ;7z-cell 
containing x. See [Str06, §2.1-§2.2] for details. 

Definition 1.8. The normal derivative of a function u is computed at a boundary point 
q e Vo by 

d^u(q) lim — y {u{q) - u{y)), q e Vq. (1.11) 

At a general junction point x - F^^q, the normal derivative is computed with respect to a 
specific m-cell K^j: 

d!^"u(x) = dl"u(F^q) := ^ 5„(moF„)(^). (1.12) 

1 .2. Statement of main result. 

Tlieorem 1.9. A^iMme f/iaf /l is not a Dirichlet eigenvalue of A, and neither is rcai-icaA, for 
any a> e W,. For the Laplacian on X with Dirichlet boundary conditions, the resolvent 
kernel G' " defined by (1.1) is given by the formula 



G'-"(x,y) = Yj '■o.'i""'-'\F-Jx, F-Jy), (1.13) 
where ^'\x,y):^ J] G%il,'^{xW^'{y). (1.14) 



p,qeVi\Vo 

where convention stipulates *i"'"''""(F^'x, F^^y) — Ofor x,y not in F^X. In formula (1.14), 

,(.1) 
p 



i^/'p is the solution to the resolvent equation at level 1, i.e. 



{Al - A)i/r'„" = 0, on each K, = F,(X), 
< J J (1.15) 

\'l^p\q) = ^/;?, for p eVi\Vo and qeVi, 

where 6pq is the Kronecker delta. The coefficients Gpg in (1.14) arise as the entries of the 
inverse of the matrix B given by 

fi» J] S^^iASf (?), qeFjm, (1.16) 

Kj3q 

where the sum is taken over all 1 -cells containing q. 
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This result appears with proof as Theorem 3. 12; a similar formula for Neumann bound- 
ary conditions appears in Theorem 4.2. 

Remark 1.10. In (1.16) and elsewhere, we use the notation YjKi3q to indicate a sum being 
taken over the set {j '-. q e Kj - Fj{X)]. 

The rationale for the definitions ( 1 . 1 3)-( 1 . 1 6) is best explained by the following heuris- 
tic argument and by comparison to [Str06, Thm. 2.6.1]. One would like to be a weak 
solution to the resolvent equation on a 1-cell C - Fi{X), except at the boundary where 
some Dirac masses may appear However, this implies that ri^"''"'"^'{F'f^ x, F'^^y) will be 
a weak solution on the 2-cell F -,{€), and in the limit (1.13) gives a solution on the entire 
fractal. Each term added to the partial sum of (1.13) corresponds to canceling the Dirac 
masses at the previous stage and introducing new ones at the next; these are wiped away in 
the limit. 

For to be a weak solution at level 1, we mean that if m 6 dom A and u vanishes on 
dX = Vo, then 



^"\x,y)(Al-^)u(y)d^i{y)^ ^ i//^\x)u{p). 

peV,\Vo 

With (1.14) as given above, integration by parts and linearity give 

r r"(x,y)(^I-A)«(y)flfM3')= r [(AI-Ay)'i"'>(x,y)]u(y)dfi(y) 

= X r [(M-A)t/^','\y)]u(y)dti(y), 

where we used the notation A, to indicate that the operator A is applied with respect to the 
variable y. 

Now by (1.15), t//q' satisfies the resolvent equation on the interior of the 1-cells, but 
-Ai//g' has Dirac masses at the boundary points with weights B^l :- YjKj3sdn^''q{s). In 
otherwords, wehaveAt^J^" - ii/ri^" except on Vi \ Vo, so that (/iI-A)i/c^,"(y) - Y^seVAVa^^'ql^siy), 
where 6s is the Dirac mass at s. Therefore, the calculation above may be continued: 

r 'V'\x,y)(Al-A)u{y)dti{y)^ ^ ilf'^^\x)G'»^B'^» f 6 Ay)u(y) dfiij) 



peVi\Vo 

The foregoing computation is the origin and motivation for (1.14)-(1.16). A key tech- 
nical point is the use of a linear combination u of vectors tfr'q' for which {AI - A)u is a 
single (weighted) Dirac mass at p. From the calculation, it is clear that this hinges on the 
invertibility of B; this is the significance of Lemma 3.7. 

As mentioned just above, once the solution is obtained on level 1, it may be transferred 
to a cell FcoiX) by rescaling appropriately. However, this is not sufficient to allow us to 
compute (AI - Ay)G''"(jc,3'); some finesse is required to ensure that these solutions match 
where these cells intersect, that is, on the boundary points Vm+i \ Vq. Some further work is 
needed; this is carried out in the technical lemmas of §3. 

2. The resolvent kernel for the unit interval 



The unit interval I — [0, 1] has a self-similar structure derived from the IFS consisting 
of Fi{x) - I and F2{x) - j + j- In this section, we exploit this perspective to derive 
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Figure 1. Mathematica plot of "F"' from Prop. 2.1 for A = I. (a) T'\x,y). (b) 
|(vj/(.i/4)(2j._23>) + 1"-"*>(2x - l,2y - 1)). (c) i(T<"""(4A:, 4v) + >P<"'«(4a- - l,4v - 1) + 
<|/(.vi6)(4-^ - 2, 4)> - 2) + 'P'"'"(4x - 3, 4)> - 3)). 



the resolvent kernel for the Dirichlet Laplacian on / by mimicking the construction of the 
Green function in [KigOl, §3.5] (see also [Str06, §2.6]). This exposition is intended to 
make the general case (presented in the next section) easier to digest. We build towards the 
result stated formally in Prop. 2.1. 

Proposition 2.1. Let A = -j^ be the Laplacian on the unit interval I — [0, 1], taken with 
Dirichlet boundary conditions. If A is not a Dirichlet eigenvalue of A, then the resolvent 
kernel G' '* in (1.1) is given by 

oo ^ 

G^\x,y) ^'i""""'(F-Jx, F-Jy), (2.1) 

H7— |tt>|=m 

sinh-^ 

for r"(x,y) ^^W(x),A*"(3'), (2.2) 

2 VI cosh 4^ 



1 I sinh VIjc, X < 



1 



and i/r'-"(x):= --2' (2.3) 

sinh-|i isinh ^|A(\-x), x > i, 

where convention stipulates *l"''''*"''(/^„'j^, P^^y) — Ofor x,y not in F^L 

Remark 2.2 (A preview of the general case). Note that the sum in (2. 1) is finite if x y, or 
if .X = y is dyadic. More importantly, i/r"* = i^j'^', is the solution to the resolvent equation at 
level 1, i.e. 

f(^I-A).A'-» = 0, on(0,i)and(i,l), 
|iA'''(0) = '/'''"(l) = 0, and lA'-'X^) = 1- 
In §3, we develop the resolvent kernel in the general case from these observations. 

In keeping with the self-similar spirit of the sequel, we use the term 1-cell in reference 
to the subintervals [0, j] and [\, 1] in the following proof. 

Proof of Prop. 2.L On the unit interval /, one has the resolvent kernel 

,,, 1 I sinh VI(1 - y) sinh Vlx x <y, 

G^-'\x,y)^— -{ ' ^' ^ (2.5) 

Visinh \A [smh yAy smh yA(l-x), x>y. 

For X < J <y, one has 

G^-(x,y) = ^i"hVI(l-y)sinhVI. ^^.5) 
VI sinh VI 
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Figure 2. Mathematica plot of G'" for A = i and two of its partial sums, (i) The sum 
of (a) and (b) in Fig. 1. (ii) The sum of (a), (b), (c) in Fig. 1. (iii) The resolvent kernel 
G''"(x, y) obtained in the limit. 



sinh sinh VIx sinh VI(1 - y) 

— Sinh 2a - 2 sinh a cosh a 



2 VI cosh ^ sinh-^ 
sinh^ 

- — ^«A'"W'A'*0') by (2.3). (2.6) 

2 VI cosh 

The same computation can be repeated for y < ^ < x and hence (2.6) holds whenever x 
and y are in different 1 -cells of /. 

It remains to consider the case when both x and y lie in the same 1-cell of /. Suppose 
that X < y < J and consider the difference 

sinh 

R(x,y) := G'"(x,y) L^^W(x)iA'''(3') 

2 VI cosh 4^ 

sinh VI.^(sinh VI(1 -y) - sinh VIj) 

VI sinh VI 
sinh VIxsinh VI(^ - y) 



(2.7) 



VI sinh ^ 

^ ^&"''(2x,2y), (2.8) 

where (2.7) follows by (2.5) and the identity sinh(l - a) - sinha = 2sinh(j - a) cosh \. 
In the case when y < x < j, one also obtains Rix,y) - jG'-'^'*'(2x, 2y). On the other hand, 
when X and y are both in the other 1-cell, one obtains (by analogous computations) that 
R{x,y) = iG'"'*'(2jc - \,2y — 1). Note that if A is not a Dirichlet eigenvalue of A, then 

neither is 1/4'" for any m - 0, 1,2, Consequently, if we define *I""(ji:,y) as in (2.2), 

then formula (2.1) for G-^Kx, y) follows. □ 

Remark 2.3. It is interesting to note that the coefficient which appears in (2.6) is 

sinh-^ 1 

2 VI cosh ^ " ^w'(i-)-<A«'(^+)' 

Formally, this indicates {Al - K)G'''^'{x,y) - 6{x - y); compare to [Str06, (2.6.3)]. Also, 
observe that 

sinh-^ 

G'\x, i) = ^ ^ ^ W(x)^'-'>(1). 
2 VI cosh ^ 



8 



lONESCU, PEARSE. ROGERS, RUAN, AND STRICHARTZ 



At each successive iteration of (2.8), one is essentially "correcting" the formula on the 
diagonal for the m-cell with rescaled copies of the formula for the (m + l)-cell; Figures 1 
and 2 are intended to explain this. In the next section, we follow this strategy for the 
construction of the resolvent kernel in the general case. 

Remark 2.4. The procedure in the proof of Proposition 2. 1 may also be carried out for the 
Neumann case: define a function i/j" " to be the solution of 

'(^I - A)^'" = 0, on [0, i] and [i, 1] 
^.<f^'>(x) = 0, X = 0, 1 



dx 



which is given by 



W'Hi) = 1 



1 (cosh y/Ax, -"^ - 5' 



cosh \cosh VI(1 - x), X > i. 
Observe that in parallel to Remark 2.3, one again has 

cosh 



2 VI sinh ■ 2 



and 



cosh-^ 



2 VI sinh f 
By analogous computations, if we define 

cosh — 

^'fcy) = ^ — ^^»(x)^<%), 

2 VI sinh 



then we obtain the Neumann resolvent kernel 

] 

2 



ni=0 \aj\=m 

3. The DiRICHLET RESOLVENT KERNEL FOR P.C.F. SELF- SIMILAR FRACTALS 

In this section, we proceed through a sequence of lemmas which will allow us to prove 
Theorem 1.9, which is stated in full in Theorem 3.12. On a first reading, the reader may 
wish to read Theorem 3.12 first, and then work through the lemmas in reverse order We 
take one hypothesis of Theorem 1.9 as a blanket assumption throughout this section: 

Assumption 3.1. None of the numbers /l^ - ji^r^A, for o) e W«, is a Dirichlet eigenvalue 
of the Laplacian. 

We construct the resolvent kernel formula according to the following rough outline: 

(1) We build a solution rfp to the eigenfunction equation which takes the value 1 at one 
boundary point of X and is on the other boundary points. 

(2) We show how i^^' may be written in terms of rescaled copies of jfp, i.e., we decompose 
the solution around a point p eV\\ Vq into solutions for each cell containing p. 

(3) We use this construction to obtain a solution on the cells of level m. 

(4) We show how the {m + l)-level solution contains Dirac masses on V,,, which cancel 
with the Dirac masses of the m-level solution, so that the sum over m is telescoping 
and yields a global solution. 
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The first two steps are carried out in §3.1. In §3.2, we collect some properties of Bp'^ :- 
YjKjBcjdn'^pXl)' as introduced in (1.15)'. For each A, we think of Bp^ as the entries of a 
matrix in p and q. Under Assumption 3.1, we show B' " is symmetric, invertible, and that 
lim^^o B' " = B"". Finally, the remaining two steps are carried out in §3.3. 

Throughout this section, we will need to analyze the properties of a continuous function 
that satisfies the /l-eigenfunction equation on all 1 -cells, but whose the Laplacian may fail 
to be in Lr. Our motivation is that the Laplacian of such a function has Dirac masses at 
points p eVi\Vo with coefficients that can be computed from the normal derivatives. The 
following result is standard; see [Str06, §2.5], for example. 

Proposition 3.2. If u is continuous and Au — vj on each l-cell Kj — Fj{X), then Am exists 
as a measure, as in (1.7), and 

J 

AM(x) = 2y^^,(x)- Yj 5,{x)Y^d''u{q) (3.1) 

where dq{x) is a Dirac mass and dn'u{q) is the normal derivative of u at q with respect to 
the cell Kj, and the sum is expressed in the notation of Remark 1.10. 

3.1. The basic building blocks of the resolvent kernel. 

Lemma 3.3. For any A that is not a Dirichlet eigenvalue of the Laplacian, and for each 
p e Vo, there is a function r['p\x) e domyvi A, as in (1.7), which solves 

{Al-A)rf^{x)^0, onX, 
ipiq) = 6pq, V<? e Vo, 

where 6pq is the Kronecker delta. Moreover, if (p is the harmonic function on X with 
^p{q) - Spq, then 

rf'^ ^i^p- AGf'^ on all ofX, and (3.3) 

dnrj'pKq) = d„(piq) - Ak';> for q € Vo, (3.4) 

where 0jf and K'pq are meromorphic functions of A with poles at the Dirichlet eigenvalues 
of the Laplacian, and Kpq — Kgp. 

Proof. Let {/„) denote the Dirichlet eigenfunctions of the Laplacian, with the correspond- 
ing eigenvalues A„ arranged so that A„+[ > A„; equality occurs iff A„ has multiplicity greater 
than one. The functions may be assumed orthonormal, and their span is dense in L^. 
Consequently we may write = 2„ ap(n)f„. The function 

_ y ap(n) 

n 

then satisfies 

(^i - A)0« = Yj jrr^'^ " " Tj "p^''^^" " 

in the L^ sense. We also see that Afljf - Y^n ^nap{n)fnl{A — A„) is L^ convergent, so 
dl'p G dovaii A, as in (1.7). In particular, is continuous and equal to zero on Vq. 
Define r[p' := - Aefp\ Then (Al - A)rf'^ = (^I - A)(p - A^p = 0, and for q e Vo, 

rfpKq) = (piq) = 5pq. (3.5) 



Recall that if/p is the solution to the resolvent equation on level 1 as defined in (1.16). 
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To verify (3.4), we will need the fact that 

which follows by computing the normal derivatives as follows: 
drJl'^iP) - d,ip\q) = 2 (/7;,"(^)5„<(^) - <'(^)5„77«(.)) 



(3.6) 



seVo 



^(77W(x)A77W(x) - i^\x)Anfix)) dfi(x) 



0, 



Now (3.4) follows via 



seVa 



Jx 



(p(xW'\x)dti(x) 



,p(x)e'g'\x) - (p{x)(,{x)) diu(x) 
^ Zj ^-^ A 2_^ap(n)aq(n) 

n ^ n 

_ ^ y, A„ap(n)ag(n) 

Define for each p,q & Vo the functions 

(.1) V '^nap(n)aq{n) 



by (3.5) 

Gauss-Green 
At?'/' = Aif. 

by (3.6) 
by (3.5) 

Gauss-Green 
by (3.3) 



pq 



Z 



A — A„ 



(3.7) 



so that d^(q{p) — d^rfqip) - AKp^. It is evident that Kpg is symmetric in p and q. It is also 
meromorphic in A with poles at the points A,,, as may be verified by writing the expansion 
on the disc of radius r centered at z (where z + An for any n, and r - inf„ |z - /l„|/2) as 
follows: 

\k 



^ ,^ _ A„ap(n)aq(n) z-A \ ^A■ V 



A„ap{n)ag{n) 



and using the fact that {«/,(;?)) and {aq{n)} are in hence their product is in f', while 
■^n/i^n - z)*^' is bounded for each k. Note that r > because the eigenvalues of a p.c.f. 
fractal have no finite accumulation point; cf. [KigOl, §4.1]. An almost identical argument 
shows that (fp is meromorphic in A with values in dom£2 A, so the proof is complete. □ 

Corollary 3.4. Let p e Vi\ Vo- Ifrj/UjA is not a Dirichlet eigenvalue for any j with p e Kj, 
then ifj and i] are related via 



\Ti^f\Ffx) ifp,xeKj 



otherwise. 



(3.8) 
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Proof. From (1.8) we have A(moF^. ') - (rjfij) ^{Au)o Fj \ for any u. Then from (1.15) 
and (3.2), one can observe that 

1/?^'" ° = ^F,(P), € Fj(Vq). 

Remark 3.5. It is helpful to compare (3.8) to the discussion of the unit interval, where (2.8) 
may be rewritten as 

[iG'"*'(2x,23;) if x,y e Kj, 
1 otherwise. 



R(x,y) = 



3.2. The matrix B*". In the construction of the resolvent kernel, the matrix Z?"' plays the 
same role as the transition matrix for the discrete Laplacian on Vi in the corresponding 
argument of Kigami for the construction of the Dirichlet Green's function. We now collect 
some important properties of " for use below. 

Lemma 3.6. The matrix B''* is symmetric for any A, and lim^^o 5*'* - B'"'- 

Proof From (3.6) we have 5n'iA^"(^) = d^'i//q\p), and thus Bj^^ = B|,'^. Then from (3.8), if 
j\, . ■ . jk are those j for which Kj contains both p and q, then 



= Z '•^'^n'f^V.) (^^'(^)) + Z by (3.4) 

1=1 (=1 -'' '' 

A- 

1=1 ■'' 

in which the final sum term is a meromorphic function of X with poles at those X for which 
rj^yLj X is a Dirichlet eigenvalue. We used the observation that the harmonic case with 
functions ^ is just the case X-Q. From (3.9) it is also clear that B^^ — » B^^ as /I — > 0. □ 

As noted in the discussion following the statement of Theorem 1.9, it is important that 
the action of B' " on the subspace Y\ \ Vo is invertible. 

Lemma 3.7. If X is not a Dirichlet eigenvalue then B' " is invertible. 

Proof. Suppose that B' " — [b^'^J is not invertible, so there are values (not all 0) 

for which 2<,ei/,\v„ B^pcflq = 0. Define 

u{x) := ^ agif/'^l\x). 

It is clear that {XI - A)m = on each 1-cell, and that =0. Now using the notation from 
Remark 1.10, we compute the sum of the normal derivatives of u over cells containing p, 
for any p eV\ \ Vq: 

Y.d'^'uip)^ J] a,Y,€'^'^\p) 

Kj3p qeVAVo Kj3p 
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3M) 

qp 

0, 



where the last equality follows by applying the symmetry established in Lemma 3.6 to the 
initial assumption. So Proposition 3.2 implies A« is continuous. It follows that {Al - A)m = 
on X, so u is a Dirichlet eigenfunction with eigenvalue A, which is a contradiction. □ 

The next result is used to prove Lemma 3.11 and also makes use of (3.8). 
Lemma 3.8. For p e Vi \ Vq and q eVowe have 

2 = (3.10) 

seV,\Vo 

Proof. For a 1-cell Kj - Fj{X), the Gauss-Green formula gives 

= r (iA5;"(x)A77»(x) - i^\x)Aifr^^\x)) dfi(x) = 
Jkj 

because both i//'pix) and rfqi^) ^re Laplacian eigenfunctions with eigenvalue A on each 
1-cell Kj. However for i e Vi we have i//''p{s) - dps, so this becomes 

dn%'\p)- 2 77«(.)5„X"(^)- (3-11) 

seFjiVo) 

The continuity of the Laplacian of 77^" at p e V\ \ Vo implies that its normal derivatives 
sum to zero, as indicated by Proposition 3.2. Thus, summing over 1-cells yields 



= Xi <'(/') = Xi Tj <(^)^n>;f(^) by (3.11) 

= ^ 'Z^" ( « ) ^ dn fp" ( s) interchange 

.seVi Kj3s 
ieVo .seVi\Vo 



where we used the sum notation of Remark 1.10. □ 

3.3. Construction of the resolvent kerneL Now that we have obtained some necessary 
properties of B' ", we can proceed with the development of a sequence of technical lemmas 
required for the proof of the main result. We begin with another corollary of Proposi- 
tion 3.2. 

Corollary 3.9. IfpeVi and A satisfies Assumption 3.1, then 

(^I-A)^w= 2 B%6,. (3.12) 

Proof. With i/r^' and B% defined as in (1.15)-(1.16), this is clear from (3.1). □ 
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Remark 3.10. From the definition in (1.16), we have = YjK,3q(^n'>J^p'(l) fovq € FjiVo). 
Thus Corollary 3.9 expresses the fact that an application of the resolvent to i/^Jf leaves 
behind nothing but a Dirac mass at every point of Vi \ Vq, each weighted by the sum of the 
normal derivatives of i//p'. 

The conclusion of the following lemma appears very technical but it expresses a straight- 
forward idea: at each stage m, our formula for the resolvent corrects Dirac masses at the 
m* level and introduces new ones at the {m + 1)*. Thus, summing over m (as we do in 
Theorem 3.12) produces a telescoping series. This makes precise the comment "these are 
wiped away in the limit" from the introductory discussion of the main result. 

Lemma 3.11. Define ^pl„, to be the unique function solving 

I (A — /i)f'„'L = 0, on all m-cells, 

I , (3.13) 
{^pAq) ^Spq^ for p e V,„ \ Vo and q e V,„. 

Then one has the identity 

(il-A,)^ r^r-"'-'\F-Jx,F-Jy) 

M=m 

Proof. Since is a sum of functions satisfying the /l-eigenfunction equation on the 

level 1 cells Kj, it is immediate that 

(A - Ay)Y^~>'~-\F-Jx, F-Jy) = 0, for y i y„,+i . 

By Proposition 3.2, we therefore need only compute the sum of normal derivatives at points 

of y„,+i. 

(1) First suppose that z e y,„+i \ V,,, with z - F^jp for some |w| = m and p eVi \ Vo, so 
that 

^^'"-'^F-jx^F-jz) = G[:r''^r'''\F-j^)K"'"'\F-jz), 

s,ieV,\Vo 

and collecting normal derivatives at z yields 
because 

^ir" - X ^"^'^'""■"(^^ by (1.16) 

Kj3p 

= 2 d'^^'r^'iF-Jz) P-F-JzeKj 

Kj3Fz'z 

= 2 5^'^''^'>r-"oF^')(z), 

F^(Kj)3z 

where the last line follows from 5n"M(F^'^,) - r^d^^iuo F^^Xq,); cf. (1.12). 



14 



lONESCU, PEARSE. ROGERS. RUAN, AND STRICHARTZ 



Continuing the computation from (3.14) and making use of G := B we have 



= r„V;;"'-'>(F^'(x)) 

thus showing that {Al — A,) has a Dirac mass ^l\^Y{x)6ziy) at z e V„,+\ \ V,„. 

(2) Next consider a point z e V,„ \Vq. In this case there are several words w, for 
which z - F^XPi) for some p, € Vq. For such a word o) and such a p we substitute from 
Lemma 3.8 into (3.14), obtaining 

Y,d:'''''\r^'v'''^-\F-jx,F-jz)) ^- Yj G!-"-v^''"^'(^'.'^)£ir^r-''(^) 

Kj3p q,s,reVi\Vo 

^- Y€r'\F-Jx)i^^;;^-Hq)- (3.15) 

qeVi\Vo 

The result is clearly a piecewise /i-eigenfunction on level (m + 1) with respect to the x 
variable, so is determined by its values on V,„+i- In each of the terms (3.15), the values 
are nonzero only at the points of Vm+i that neighbor z in F^(X), and they are easily seen 
to coincide with - ^['i, at these points. Summing over all cells, we conclude that at 

each z e V,„ \ Vo the operator {Al - A) has a Dirac mass - ^z',m) ^ziy)^ and the result 

follows. □ 



(3.16) 



Theorem 3.12. Let if/p' be the solution to the resolvent equation at level 1, i.e. 

{Al - A)t//'p' = 0, on each Kj = Fj{X), 

>f^p\q) = Spg, for peVi\Vo and qeVi, 

where 6pq is the Kronecker delta. 
Define the kernel 

&\x,y) = Yj rJi'"*-\F-Jx,F-Jy), (3.17) 
where 'i>'\x,y):^ Y G%>lf';,\xW^'{y). (3.18) 
The coefficients Gp^ in (3.18) are the entries of the inverse of the matrix B given by 

Kj3q 

the sum taken over all 1 -cells containing q. 

For A satisfying Assumption 3.1, G*'"(.x,3') is symmetric and continuous in x andy, and 
is in domx A,, with {Al — A^,)G''"(.*^, y) — Sx{y). As it vanishes on Vo, it is the Dirichlet 
resolvent of the Laplacian. 

Proof. The symmetry of G''^\x,y) is obvious. Next, note that 

M 

(il-A,)^ Y rJV^'"'^\F-Jx,F-Jy) 

m=0 \tij\=m 
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in=0\peV„+,\VB peV,„\Vo 



by Lemma 3.11, so that 



lim (^I - A,) y y r^»P<-''""(F„'^,F^'y) = ^.vW, 

H7=0 |(i»|=n; 

in the sense of weak-* convergence. It follows that G^"(jc,y) is in domAt(Ay) and that 

All that remains is to see that G''"(j^,3') is continuous. However, Lemma 3.6 shows 
Bpg"''^ B'pg as |w| oo, and hence Gpg'^' Cfq. In a similar manner, the relation 
T/p"''"'" - + r^fii^Adfp"''"^^ from Lemma 3.3 shows that 7^^"''-" — > as |c<j| ^ 00; in 
particular we find that 1/'^"''"'" i/f™, and the latter is piecewise harmonic and bounded by 
1. The conclusion is that ^"''"|'''■^' is bounded as — > 00, and since r^j is a product of \a)\ 
terms, all of which are bounded by max, r,- < 1, 

00 

H7-0 |(it|=/j; 

is bounded by a convergent geometric series. Note that, for each m, only a finite number 
of terms in the second sum are nonzero. As all terms are continuous, so is G''". □ 



4. The Neumann resolvent kernel for p.c.r self-seviilar fractals 
In Theorem 4.2, we give the formula for the Neumann resolvent kernel. 
Lemma 4.1. If A is not a Neumann eigenvalue then there is Cpl such that 



2 C«c?„/7<-"(x) = 6p 

qeVo 



for X e Vo, <^nd Cp^ is symmetric in p and q. 

Proof. Since A is not a Neumann eigenvalue, the set of vectors |(i9„;7!y'"(jii:)) | is lin- 
early independent, whence the existence of the Cp^ is immediate. Symmetry follows from 
(3.6) because the matrix [cj^'^j is the inverse of the symmetric matrix |^c)„77p"(^)j. □ 

From this and Theorem 3. 12 we may readily deduce the following result. 
Theorem 4.2. If A satisfies Assumption 3.1 and also is not a Neumann eigenvalue, then 

&»{x,y) = &\x,y) + 2 C^^(x)Tf^(y) (4.1) 

is symmetric, is in domyv((Ay), and satisfies {A — hy)G''^(x,y) — 5Ay^ on X \ Vq. It has 
vanishing normal derivatives on Vo (^nd is therefore the Neumann resolvent kernel of the 
Laplacian. 
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Proof. The symmetry of G^'^{x,y) is immediate from the symmetry of G^'^\x,y) and of 
C%. Both G'\x,y) and if-^iy) are in domA((Ay) and {A - AyWp'iy) = on X \ Vq so 
(A - Ay)G'^\x, y)^{A- A,)G'\x, y) = 6 Ay) on X\Vo. 

It remains to prove the assertion about the normal derivatives. We will use the notation 
(5n)yG"' for the normal derivative of G'-'^\x,y) with respect to its second variable. Since 
G*'''(x,y) e domyviCA,.) it has a normal derivative at p e Vo, and by the Gauss-Green for- 
mula, 

{dn)yG''\x,p) = 2 {(dn)yG'\x, S^^^s) - G'\x, s)dnll)'M) 

seVo 

((A,G«(x, s)) 77«(.) - G'\x, s) (A,77«(i))) dp(s) 

= r (A,,-A)G''\x,s)Tjf{s)dfi{s) 
Jx 

= (4.2) 

where at the first step we used that G'"(x, s) = for i e Vo and at the last step we used 
(Aj - A) G''^\x, s) = -Sx(s) as a measure. It follows that at each p e Vo, the normal 
derivative of (4.1) vanishes: 

(d,)yG';^'(x, p) = -i'^{x) + (5„)y 2 0^li^\x)ijf{p) by (4.2) 

= -rf^Xx) + 5q_,if^\x) by Lemma 4. 1 

?eVo 

= 0. □ 




5. Example: the StERPiNSKi gasket SG 



Recall the harmonic extension algorithm as described in [Str06, §1.3]: if the values of 
a function u are specified at the points of Vo and written as a vector 



I Vn 



u(po) 
u(pi) 

uipi) 



then the harmonic extension of u to Fi(Vo) (the boundary points of the 1-cell Fi{SG)) is 
given by 



: AlU. 



Vo 



u{Fipo) 
u(Fipi) 
uiFiPi) 



where 



5 










1 

" 5 


" 2 


2 


1 


2 


2 


1 







5 





2 


1 


2 




1 


2 


2 



and At = - 



2 1 
1 2 




are the harmonic extension matrices. In general, mL „ = A„m ,, , where = A^ • • • A^j, . 
Thus, the harmonic extension matrices allow one to construct a harmonic function with 
specified boundary values. Similarly, spectral decimation provides matrices which allow 
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one to construct an eigenfunction with specified boundary values. For example, 

, [ (5-A){2-A) 

Ao(A) = 



(5 - A)(2 - A) 



(A-X) {A -A) 2 (5.1) 
(4 - ^) 2 (4 - A) 

is the analogue of Aq - Ao(0). By the usual caveats of spectral decimation, these extension 
matrices can only be used when A is not a (Dirichlet) eigenvalue. 

Remark 5.1 (Spectral decimation). A very brief outline of the method of spectral decima- 
tion is as follows. 

(1) Begin on some level m - mo with u„, and A,,, that satisfy -A,„m„, = /1,„m„, on V„, \ Vq. 

(2) Extend u,„ to a function u„, on V„,+i \ Vq by comparing the eigenvalue equations from 
each level. 

(3) Obtain a collection of extension matrices, one for each mapping in the original IFS, 
and a rational function g which relates the eigenvalues on one level to the eigenvalues 
on the previous level by g{A„j) = /l,,,-! . 

(4) Inductively construct a sequence {A,„^^, A,„^+i, A„^+2, . • . ) by choosing A,„+i from the set 
g^\A,„) for each m. 

For every such sequence that converges, alim„^cx,P"'A,„ will be an eigenvalue of A on 
X, where a and /J are constants specific to X. For the Sierpinski Gasket SG, a — \ and 
yS = 3| =5. Note that the calculations in (2)-(3) will forbid certain choices, so some care 
must be taken in the construction of {/!,„). See [Str06, §3] for more details. 

To obtain the numbers Z?J;'^ (appearing in (1.16)) for the Sierpinski Gasket SG, we find 
the normal derivatives of the eigenfunction that has boundary values (1,0, 0), as computed 
at each point of Vq. If (-il - A)m = but u is not a Dirichlet eigenfunction, then consider u 
on Fq (Vo). By spectral decimation, this is given by 

where the matrix Aq{A) is as in (5.1). We actually only need the values of the normal 
derivative 

d^uiqi) = lim (^) {2u{qi) - m(F;"^,_i) - uiFfq^)) , qi e Vq. (5.2) 

The factor | arises here because ry - | for each j in (1.3) for SG\ see also (1.12). The 
calculation of rj = | is given in [Str06, §1.3]. 

It is extremely easy to compute the normal derivatives of a harmonic function: one 
does not need to compute the limit, as all terms of the sequence are equal; see [Str06, 
(2.3.9)]. Therefore, our approach is to obtain a harmonic function which coincides with 
u on F™(Vo). The limit of the normal derivatives of these harmonic functions will be the 
normal derivative of u. An alternative interpretation would be to interpret the harmonic 
functions on SG as the analogue of the linear functions on /. Consequently, the tangent to 
a point of SG should be given by a harmonic function plus a constant, provided the tangent 
exists. This is the motivating idea of [DRS09]. 

Multiplication by Aq™ allows one to find the required harmonic function at stage m; 
rewriting the normal derivative (5.2) in vector notation, one has 

0" (2, -1,-1) (2, -l,-l)-Ao'"«|,„,,„, 

= (2, -1,-1) ■ Ao'«Ao(^,) ■ ■ ■ Ao(^i)m|^,^ . 
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It therefore suffices to understand the Umit limmA^"'Ao(A„,) ■ ■ ■ Ao(/li); this was computed 
in [DRS09]. The following theorem is the main result of [DRS09], taken with mo = 0. 

Theorem 5.2. Let a = (0, 1, 1)^, /3 = (0, 1,-1)^, y,n = (4,4 - /l„„4 - A„y. If neither of 
the values 2 or 5 occur in the sequence A^, then 



lim Aq* ■ Ao{Ao+k) ■ ■ ■ Ao(^o+i) « = T wo i n 



4A 



n('-¥)« 



limA/-Ao(.lo+i)---Ao(/lo+i)yS 

<r->oo 3 ■ J /to 



2A 



lim Aq* ■ Ao(.io+A.) ■ ■ ■ Ao(/to+i) ro = (4, 4, 



In particular, this can be used to get the desired normal derivative. We know that all we 
need do is compute 



(2, -1, -1) ■ (limAo"'Ao(.iJ ■ ■ ■Ao{A,)\u\ . 

\ in J ^0 



(5.3) 



The boundary data is be taken to be (1, 0, 0) when computing the normal derivative at 
the point p where u{p) - 1, and (0, 1, 0) at a point where u(p) - (these two points are the 
same by symmetry). Writing 



1 



we find that 



4 

A-Ao 
4-io 



A-Ao 



and 






1 





1 





1 




1 




1 





~ 2 


1 


^2 


-1 



limAo'"Ao(/l„,)---Ao(/ii) 
A-Aq 



AA 



A 3io(2 - A,) 



and by (5.3), the normal derivative is 

2(4 - Aa)A 
^Aq{2-A,) 



The normal derivative at the other point is computed by first finding 

C 

UmAo'"Ao(.l„) • ■ ■Ao(.li) 1 



1 



AA 



2 3Ao{2 - Ai) 



n('4) 



1 2A 
^2JM) 




1 

-1 



and then taking the inner product with (2,-1,-1), which cancels the second vector to leave 
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Figure 3. The 1 -cells of SG3. 



It seems logical at this point to define a function 

4-A 



riA) 



3Ao{2-Ai) 



n 



(5.4) 



and to write the normal derivative at the point where the 1 occurs as (4 - A())t(A)/2 and that 
at the point where the occurs as -t(A). We note that for a non-Dirichlet eigenfunction, 
none of the values 2, 5, 6 occur in A„, for m > 1 so the term (2 - Ai) in the denominator 
cannot be zero. It follows that 3 does not occur for m > 2 and therefore that t(/1) in 
this case. An exception to our formula as currently written occurs when A - 0, because 
then also Aq - 0, but the function t{A) is easily shown to have a continuous extension to 
A - Q with t(0) = 1; cf. [DRS09]. With this correction, our formula is also valid for the 
harmonic case. 

It is now easy to write the entries of the matrix appearing in (1.16). The term Bp], 
has two copies of the normal derivative (4 - Aq)t{A)I2, and the term B^'p^ has a single copy 
of -t(/1) at each q that is not equal to p. Both are on 1 -cells rather than the whole of 
SG, so there is an extra factor 5/3 in their normal derivatives. As a result, the matrix is 



B ■■ 



(4 - Aa)T{A) 
-T{A) 

-T(A) 



-T{A) 
(4 - ^o)t(^) 

-T{A) 



-T{A) 
-T{A) 

(4-/lo)T(/l) 



and we should invert this to get the matrix Gpq for the Green's function. Since 



det 



abb 
b a b 
b b a 



= (fl - bfia + 2b), 



the matrix B is invertible iff /Iq 5^ 2, 5, in which case 

' (3-/lo) 
1 
1 



5(5 - Aq){2 - MMX) 



1 1 

(3 - ^0) 1 
1 (3 - .lo) 



Note that this is consistent with the harmonic case where /Iq = and t(0) = 1 gives factors 
9/50 for G% and 3/50 for G% with p + q; see [Str06, (2.6.25)]. 

6. Example: SG^, a variant of the Sierpinski gasket 

6.1. The Laplacian on SGt,. The fractal SG^, is obtained from an IFS consisting of 6 
contraction mappings, each with scaUng ratio i, as indicated in Figure 3. The details of 
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Figure 4. The eigenfunction extension on one (m - l)-cell of SGj, to 
m-cells. The values on (m - l)-cell are 1,0, 0. 



the spectral decimation method for SG3 have been worked out independently in [BCD^OS, 
DS07, Zho09a, Zho09b]. Note that po is contained in three 1 -cells of SG^,, in contrast to 
each of the other points /?, of Vi \ Vo, which are contained in two. For this reason, we define 
the graph Laplacian on SG3 as 

A„Mx) = — ^ y {u(y) - u{x)), (6.1) 
deg(x) 

where deg(x) is the number of ;7z-cells containing x. From [Str06, §4.4], we have 

A^u(x) ^ lira r'" hfUfi'j deg{x)A,„u(x). (6.2) 

The renormalization constant r will be computed in §6.3. 

Let P f denote a vertex where u takes the value x as depicted in Figure 4, then by (6.1), 
the symmetric eigenvalue equations on V,„ are 

KiUiPx) ^ \ - x) + {x - x) + {w - x) + (y - x)] ^ -A',„x 

A,nu(py) ^ \ [{x - y) + (w - y) + {z - y) + {0 - y)] = -A'„j 

Amuip,) ^\[(y-z) + (w-z) + (z-z) + (0- z)] = -A'„^z 



A,„m(/?„,) = i [2(x -w) + 2(y-w) + 2(z - w)] 



(4 - A„,)y - X + z + w 



which can be rewritten, using A,„ - 4/1,',, as 

(4 - Am)x - I + X +y + w 

(4 - A„,)z ^y + z + w (4- A„,)w ^^(x+y + z). 

For now, we suppress the dependence on m for convenience and denote A = A,„. Solving 
for A, we obtain 

(6.3a) 
(6.3b) 
(6.3c) 
(6.3d) 
(6.3e) 



X = 


a(A) 


= (96- 


109A + 33A^ -3A^)/ip(A), 


y = 


m 


= (16- 


3A)(3 - A)/^{A), 


z = 


y(A) 


= (36- 


lA)/<p(A), 


w - 


pW 


= 4(5- 


A)i3 - A)/<p{A), 


where 


ip(A) 


= 3(5- 


A)(3 - A){4 -6A + A^), 



and we see that the forbidden eigenvalues are 3, 5, 3 + V5. 
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Figure 5. The labeling for a general eigenf unction extension, from one 
(m- l)-cell to six m-cells. The values on the boundary of the (m- l)-cell 
are a, b, c. 





l^{u{x)-u{y)) 

Figure 6. The harmonic extension of u on SGt,, where m|v„ = [1,0, 0]. 



For a general function on SGi, we extend the eigenfunction using the labeling indicated 
in Figure 5, as follows: 



x' = aa(A) + cl3(A) + by{A) 
y = ba{A) + ap{A) + cy(A) 
z' = ca{A) + bpiA) + ay{A) 



X = aa(A) + b/3(A) + cy{A) x' = aa{A) + c/3(A) + by(A) (6.4) 

y = ba{A) + c/3{A) + ay{A) 
z = ca{A) + a/3(A) + by(A) 
w — (a + b + c)p{A). 

The eigenfunction extension matrix for SG3 corresponding to Fq is 

1 
Aq(A)^ aiA) m y(A) 
aiA) yiA) p{A) 

where we have a{A),l3{A),y{A), (p(A) as before, that is Ao(A)u\v„ = mIfoVo- 

6.2. Eigenfunctions of the Laplacian on SG^. Let m be a function taking values 1,0,0 
on Vq. The harmonic extension u on V] corresponds to taking /I = in the system (6.3) 
above, so that 



15' 



y = 



15' 



15 



w - 



15' 
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1 

1 ~ f- — \ 1 ~ ^m+\ 




Figure 7. The values of the eigenfunction u on ^Gs, where u\\/^ - 
[1,0, 0]^. This figure shows a closeup of u near the point where it takes 
the value 1. By symmetry, we define x,„ := 1 - u(F'^q\) = 1 - u{F'l^qi). 

and we have Figure 6. Following [Str06, §1.3], the energy renormalization constant com- 
puted by 

= (4 ■ 1^ + 2 ■ 2^ 4 ■ 3^ 4 ■ 4' + 2 ■ 7") = 

Since E(,{ii) =1 + 1=2, 

2 = £o(M) = r->£i(fi) = r~i|^ =^ r=^. (6.5) 

Thus, the normal derivatives on SG^, are computed by 

/ 15\'" 

5„m(p) = lim — y {u(p) - u{fi). (6.6) 

p~my 

Theorem 6.1. The pointwise formulation of the Laplacian on SG^ is 

/90V" 

Af,u(x) = 6 hm I — I A,„u(x„,), (6.7) 



in—^oo 



, 7 , 

where {x„,} is any sequence with limx,„ = x and x^ e V, 
Proof. Following [Str06, §2.2], it is easy to compute 



if deg(x„,) = 4, 



g,„ ;/ deg(x,„) = 6 

since fi is the standard (self-similar) measure on SG3. Thus, by (6.2), 

/ISV" /90\"' 
A^u(x) = Um — ■ 6"'+'A,„m(x) = 6 lim — A,„m(x„,), (6.8) 

m— >oo \ I j m^co \ I j 



Throughout, whenever there is discussion of an eigenvalue A, we assume that we have 
been given the sequence {/l,„))^^o which defines A via the decimation formula. Thus by 
Theorem 6.1, 

/90r , 3 /90\"' 
i = 6 1im - ^; = -lim - A,„. (6.9) 

m— *oo \ I j L w?— >oo \ / I 
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6.3. Computation of the normal derivatives. 



Theorem 6.2. Let -Am = Au on SG^, where u is defined on Vq by uiqo) — 1, — 0, 

and u(q2) — 0. Define 



... lift 



3^0 1 



2^J + T 



Then the normal derivatives of u are 



d„uiqo) = — ^-^t(/1), and 
d„u{qi) = d„u(q2) = -t(.1). 



(6.10) 



(6.11a) 
(6.11b) 



Proof. To obtain (6.11a) we need the values uiF^qy) - uiF'^qi) - \ - x,„ as depicted in 
Figure 7. 

Claim; u{F'^'qi) = u(F'^q2) - \ - x„, where xq - 1 and 

^m+l (4 - /tm+l)(6 -/!,„+ 1 )/!„,+ 1 / A„ 



(4-6i„,+i +i2^j)^ 



(6.12) 



Proof of claim. By (6.4), if = c, then x - x' . Then from u{qi) - uiqi) = 0, we will 
have uiFg'qi) = u{F'^q2) for all m, by induction. Define x„, := 1 - uiF^'qi), m > 0. From 
u(qi) - 0, we have xo = 1. Now we show (6.12) holds. 

Denote 6(A) /3{A) + y{A), so that 6{A) = (14 - 3A)(6 - A)/(p{A), where (p{A) is as in 
(6.3). Using (6.4), we have the matrix equation 



gives 





1 




1 


AoiA) 


1 ~ -^m 




1 - x,„+i 




1 ~ -^m 




1 - X,„+l 


Xm+l = 1 - 0!(A,n+i 


) - ^(^i+I 


) + S(A,„+i)x,„ 



^m+\(5 - Am+l) _^ (14 - 3/l„,+ i)(6 - A„+i) 



4 - 6/l„+i + /l^ 



From the decimation relation [DS07, (2.12)], we have the identity 

3(5 - /lm+i)(4 - A,„+i)(3 - /l„+i)/l„+i 



so that 



SiAm+i) = S{Am+i) 



{U-3A,„^i)A,„ ^' 

3(5 - /l„,+ i)(4 - /l,„+i)(3 - Am+l)Am+l 



{U-3A„^i)A„ 

_ (4 - A,„+i){6 - /tm+Q/^m+I 

(4-6/l,„+i +^2^j)/},„ 

We would like to see x,„+i - /(/Im+i) = ^(/lm+i)(j^m - fiAm)) for some function /, which is 
equivalent to 

(4 - 6^,+i + ^2^1):^*^^ = Z^i:;!2(4 _ ^,,,^;)(6 - A,„^,) - (5 - ^,„+,). 

1 

Let f{x) - xg{x) and this can be rewritten 

(4 - 6^,+i + Al^^)g{A,„^i) = g(A,„){24- - 10^,+i + /l^^i) - (5 - /l,„+i). 
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which is easily seen to be true for the constant function g(x) - | . Hence we may define 
f{x) = |, to obtain 



Now we compute d^u{qo) using (6.12) to obtain 



4-Ao 



4^0 



A) 



(4-Aj)(6-Aj) 
4 - 6Ai + A^. 

■' J 



{A-Aj){6-Aj) 
4-/l() ' i.^f A-(,Aj + A] 



n 



Since m(^o) = 1, we apply (6.6) to compute 



15 



4-^0 
2Ao 



d„u(qo) = lini^ ( — ) I - I (2u{qo) -2(1- x,„)) 



lim I A„ 



Ao 



(4-Aj)(6-Aj) 
6"'(4 - Aq) ■ I f 6(4 - 6Ai + Ah 



n 



'-^I'-A 

2Ao - 



which is equivalent to the result. 



(4 - Aj){6 - Aj) 

y 6(4-6^; + ^])] 



Now we compute the normal derivatives (6. 11a). To obtain d^^u{qi) - d^u{q2), we don't 
actually need the values u{F'"q()) and u{F"'q2) as depicted in Figure 8. Instead, it suffices 
to only compute their sum, since by (6.6), one has 



/15\ 

d„u(qi) = d„u(q2) = - ]m y 1 {uiF'l'qo) + uiF'J'qo)) . 

To exploit this symmetry accordingly, define 

y,„ := uiF^qi), Zm := uiF^qo), and i„, i^y,,, +z,„. 
Claim: the sequence {sm]°^^^Q is given recurrently by sq = 1 and 

_ (14-3^,+i)(6-^,+i) 

1 — / 1 \ 

'P{Am+l) 



(6.13) 



(6.14) 



Proof of claim. As indicated in Figure 8, dihedral symmetry allows us to continue using 
the same matrix Ao{A) for computations, as long as we use [0, 1,0]^ for the new boundary 
data. 

It is clear that io = 1 + from the values on Vq. Then using the notation d(A) - 
a{A) + /3(A) as above, the matrix equation 




















AoiA) 






/i{Am+l)ym 


+ 7(Am+l)Zm 




Jin+l 








_ 7(Am+l)ym 


+ /3(A,n+l)Zm 




Zm+ 1 



gives = y,„+i + Zm+i = 6(A,n+i)Sm immediately. 
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'771+ 1 "^777. 



Figure 8. The values of the eigenfunction u on SGt,, where u\vo - 
[0, 1,0]^. This figure shows a closeup of u near a point where it takes 
the value 0. See (6.13) and the ensuing discussion. 



Since (6.14) gives 



A) ^. I 4- 6A: + A] 

and uiqi) = 0, the normal derivative is 

15V'76>"' 



d^uiqi) = hjn (2m(?i) - s,„) 

1 ,. /90V" ^ A (4-Aj){6-Aj) 



(90\ 
lim — A,„ Pf 

0)— >oo y / y y ^ 



Ao \1 I ""' A. f 6(4 - 6/1/ + A^.) 

y-i ■' J 

1 /2_^\pj (4-^,)(6-i,) 

j\j 6i4-6Aj + Aj)' 



6.4. The resolvent prekernel. As in (1.16), let fi^'^ d^'i//f\q) for /? e Vi \ Vo- 

Corollary 6.3. With t(A) as in Thm. 6.2 and r - j^, 



B« - -r-\{A\ forp~,q and B« = -> ^ 



/'P ^-i 



\-\4-Aq)t{A\ p^po 



\4-Ao)t(A), p^po. 

Proof. We are now working on Vi, so each term has a leading factor of r"'. Whenever 
p ~i q, there is just one term d^u(q) - -t{A) in the sum; the other corner of the triangle is 
ignored and everything outside this 1-cell is 0. When p - q, then there is a sum of terms 
d^u{p) - ^^t(A). At the center point po, there are three such terms; at every other point 
there are only two. □ 
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The matrix B^'pg is 



1(4 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


(4 - .io) 


-1 











-1 


-1 


-1 


(4 - ^o) 


-1 











-1 





-1 


(4 - .!,)) 


-1 








-1 








-1 (4 


- '!()) 


-1 





-1 











-1 (4 


- 


-1 


-1 


-1 











-1 


(4-/! 


Define the resolvent prekeniel by G' " 


:= (B<'") 


-1 





Our final result may be obtained by brutal and direct computation. 
Theorem 6.5. The resolvent prekeniel G' " is given by 



14 



15(6 - A)T{A)>f{A) 



(2 - A)Ki 


Kl 


Kl 


K\ 


Kl 


K\ 


Kl 


Kl 


Kl 


K3 


K4 


K5 


K4 


K3 


K\ 


K3 


Kl 


K3 


K4 


K5 


K4 


K\ 


K4 


K3 


Kl 


K3 


K4 


K5 


K\ 


K5 


K4 


K3 


Kl 


K3 


K4 


Kl 


K4 


K5 


K4 


K3 


Kl 


K3 


Kl 


K3 


KA 


K5 


K4 


K3 


Kl . 



where 



Kl 


= (3- 


^)(5 - A){6 - A), 


Kl 


= 201 


-3Q0A+'=fA^-24A' + lA' 


K3 


= 87 - 


75A + \9A^ - 


K4 


= 57- 


2AA + and 


K5 


= 51 - 


\5A-A^. 



In particular, G* '* is symmetric and invertible with determinant 

^7 ... . -,2^ 



detG'" = — 



7 \' 6(4-6A + A^) 



15/ (6 - A)ip{A)^T{Ay 
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